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Abstract 

We investigate the properties of hybrid stars consisting of quark matter in 
the core and hadron matter in outer region. The hadronic and quark matter 

x: 

H , equations of state are calculated by using nonlinear Walecka model and chiral 

colour dielectric (CCD) model respectively. We find that the phase transition 
from hadron to quark matter is possible in a narrow range of the parameters of 
nonlinear Walecka and CCD models. The transition is strong or weak first order 
depending on the parameters used. The EOS thus obtained, is used to study the 
properties of hybrid stars. We find that the calculated hybrid star properties are 
similar to those of pure neutron stars. 



1 Introduction 

The first model of neutron stars was proposed by Tolman, Oppenheimer and Volkoff 
[|IJ] in late 30 's. In this model, the properties of neutron stars were calculated by assuming 
that the neutron matter consists of non-interacting degenerate gas of neutrons. Since 
then, a number of improved calculations JJ|~@? which use different nuclear equation 
of state (EOS), have been performed to determine the composition and properties of 
neutron stars. The subject of neutron stars received a new impetus after the discovery 
of pulsars ]7|, |S|. In spite of two decades of work, the EOS required in the calculation of 
neutron star properties is uncertain for a number of reasons. The EOS at large baryon 
densities sensitively depends on short range nucleon - nucleon interaction and this is one 
source of uncertainties. Further more, at large baryon densities, weak decay of neutrons 
into hyperons may be energetically allowed and the neutron matter at high densities 
may contain hyperons. Efforts have been made to include hyperon degrees of freedom in 
nuclear EOS |?J-[I2|. However, it must be noted that the interaction of hyperons with 



other baryons is not well understood and the EOS is, therefore uncertain. 

The quark structure of hadrons implies that at sufficiently large nuclear densities 
the nuclear matter should convert itself into quark matter. The question of transition 
density and the order of the transition of the nuclear matter to quark matter is not yet 
settled. The density at which the transition occurs, is believed to be few times nuclear 
matter density. The lattice calculations indicate that for nonzero quark masses the phase 



transition may be weak first order or second order | I5| . Most of the model calculations 
assume it to be first order. Thus, for large enough mass of neutron star, its core may 
consist of quark matter. In addition, if the phase transition is first order, a part of 
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the core may consist of a mixed phase of quark and nuclear matter. Particularly, if, at 
the phase transition, the discontinuity in the baryon density between quark and baryon 
matters is large, a substantial fraction of neutron star may be in the mixed phase. It is 
therefore, interesting to investigate the properties of neutron stars having a mixture of 
quark and nuclear matter. 

In the present work, we have studied the quark-hadron phase transition at zero 
temperature using two phase model and applied it to the calculation of neutron star 
properties. Ideally, one would like to investigate this phase transition using same model 
in both the phases. One such method would be the lattice QCD calculations. However, 
the problem of inclusion of dynamical fermions in lattice calculations is not solved and 
the calculations of equation of state for nonzero baryon chemical potential will not be 
available in near future. Thus, a number of calculations employing different models in 



two phases have been done 11 0, |lQf. Here, we want to present a calculation where 



we have used nonlinear Walecka model ]TD|, [TT|, for nuclear matter and chiral colour 



dielectric model(CCDM) for quark matter ]Tjj] . 

The nonlinear Walecka model has been used extensively in nuclear structure calcula- 
tions. It has been observed that, with nonlinear self-interaction of the a field, one obtains 
very good agreement with the properties of nuclei over a wide range of the periodic table 
[ |T6|j . Thus, it seems reasonable to use such a model to describe the hadronic phase in a 



region where nuclear densities are not too large in comparison with the nuclear matter 
density. The equation of state calculation in hadronic phase includes hyperons as well as 
leptons. The parameters of the nonlinear Walecka model are determined by fitting the 
properties of nuclear matter. However, the coupling of a, p and u mesons to hyperons 
is not determined by this procedure. We have therefore varied these couplings within 



a reasonable range and investigated the effect of these on the phase transition and the 
properties of stars. 



The CCDM has been used earlier in baryon spectroscopy [|17|. These calculations 
have shown that the model is able to explain the static properties of light baryons 
very well. Further more, when applied to quark matter calculation, the model yields an 
equation of state which is quite similar to the one obtained in lattice calculations for zero 
baryon chemical potential. In particular, it shows that the energy density calculated in 
CCDM is close to the energy density of free quarks and gluons and the pressure decreases 



rapidly when the temperature is close to the transition temperature [I3|. We therefore, 
think that CCDM is a qualitative improvement over the bag model, which is often used 
to calculate the equation of state in quark matter. The model has also been used in the 
calculation of quark star properties [[L8| . 

The results of our calculations can be briefly summarised as follows. We find that 
the hadron to quark matter phase transition occurs for a narrow range of the parameter 
sets of the models of the two phases. Particularly, the phase transition does not occur 
for smaller values of nuclear compressibility as well as for certain parameter sets of the 
CCDM (see later). This probably indicates that certain range of parameters of nonlinear 
Walecka model and/or CCDM are physically not acceptable. The neutron stars in our 
model, called hybrid star, may consist of quark core, a mixed phase region and then outer 
part made up of neutrons. The width of the mixed phase depends on the parameter set 
used. Usually a stiffer EOS gives higher mass or radius for a pure neutron star. But for 
a hybrid star, trend is opposite. The neutron stars in our model is found to be consistent 
with the observational limits. 

The paper is organised as follows. In Section 2, we obtain the equations of state for 



nuclear and quark matter. Section 3 is devoted to the discussion on hadron-quark phase 
transition. The results for the neutron, quark and hybrid stars structure are discussed 
in Section 4 and the conclusions are presented in section 5. 

2 Equations of State 

2.1 Hadrons 

The equation of state for hadrons is calculated in the frame work of mean field theory 



using the nonlinear Walecka Lagrangian given below |K| [12] 
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The Lagrangian in eq(l) above includes nucleons, A and S~ hyperons (denoted by 
subscript i), electrons and muons (denoted by subscript I) and a, u and p mesons 
(given by a, 10^ and p a,At respectively). The Lagrangian includes cubic and quartic 
self-interactions of the o field. The meson fields interact with baryons through linear 
coupling and the coupling constants are different for nonstrange and strange baryons. 
The Lagrangian defined above (without strange baryons) has been used extensively in 
nuclear structure calculations with success. The parameters of the model are meson- 
baryon coupling constants, meson masses and the coefficients of the cubic and quartic 
self- interactions of a meson ( b and c respectively). The u and p masses are chosen to 
be their physical masses. 



The equation of state is obtained by adopting mean field ansatz. Thus in presence 
of baryons, the mesons develop nonzero vacuum expectation values (er, uo and p a re- 
spectively). Assuming that the baryon densities are uniform, one finds that the time 
components of uj and p 3 , in addition to a, are nonzero. One can then define effective 
masses (fhi) and chemical potentials (fa) for the baryons as, 

rhi = mi- g ai o (2) 

and 

fa = Hi- gw& - hg P NP 3 (3) 

where I3 is the value of the z-component of the isospin of baryon i. The Fermi momenta 



(fcj) and number densities (rii) of the baryons are given by k~ = J/j,? — mf and n. 
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For leptons, the Fermi momenta and number densities are given by k\ = J /jif — ml 
and m = A. 

' ■ATT'' 

The condition of equilibrium under weak interactions (assuming that neutrinos are 
not degenerate) give the following relations between baryon and lepton chemical poten- 
tials 

r"P r"n r^ e 1 r^A H"ni 

A*E- = Hn + Hp, Hfi = He (4) 

and charge neutrality gives, 

n p = n e + n^ + n s (5) 

In addition, ns = n n + n p + n^ + ns is the total baryon density and jib = Hn is defined 
as the baryon chemical potential. 



The mean field values of a, u>o and pjj are determined by minimizing the energy at 
fixed baryon density. Then the EOS is calculated using the expression for pressure P 
and energy density E as given below. 

p = 2 m ^o + 2 m ^o) 2 - 2 m ^ 2 ~~ 3 &rr M#^) 3 - l c (9«Nv) 4 

+ J2 P FG(rni,fii) + J2 P FG(mi,Hi) (6) 
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E = 2 m ^o + 2 m ?(^o) 2 + 2 m ^ 2 + ^ bmN ^N^f + -c{g aN af 

+ J2 E FG(rni,fXi) + ^E FG (m h [n) (7) 
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In the above equations Pfg an d EpQ are the relativistic non-interacting pressure and 
energy density of the fermions. Nonlinear Walecka model has eight parameters out 
of which five are determined by the properties of nuclear matter. These are nucleon 
coupling to scalar (— ), vector (—) and isovector mesons (— ) and the two coefficients 
in scalar self interaction i.e. b and c . These are obtained by fitting saturation values 
of the binding energy/nucleon (— 16MeV ), baryon density (0.15/m -3 ), symmetry 
energy coefficient (32.5MeV ), Landau mass (0.83mAr). The nuclear compressibility is 
somewhat uncertain and therefore, we have varied it between 250MeV — 350MeV. The 
values of these parameters along with compressibility are presented in Table 1. 

The other three parameters are coupling constants of hyperon-meson interaction and 
are not well known. These cannot be determined from nuclear matter properties since 
the nuclear matter does not contain hyperons. Further more, properties of hypernuclei 
do not fix these parameters in a unique way. In the literature, a number of choices have 
been made. These along with our choice are listed below. 



1. choose them to be same as nucleoli- meson coupling constants (Universal coupling 
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2. choose them to be a/2/3 times nucleon-meson coupling constants |]TTJ, 

3. choose them to be 1/3 times nucleon-meson coupling constants [I2| 



4. It is claimed that the second choice follows from quark counting rule. However, 
using SU(6) quark wavefunctions of baryons and assuming that the mesons do 
not couple to strange quarks, we find that the ratio of meson-hyperon and meson- 
nucleon couplings are 2/3, 2/3 and 1 for a, uo and p mesons respectively. 

In our calculation, we have used the above couplings to study the nuclear EOS and 
star characteristics. 

2.2 Quarks 



The Colour Dielectric Model (CDM) is based on the idea of Nelson and Patkos [19]. In 
this model, one generates the confinement of quarks and gluons dynamically through the 
interaction of these fields with scalar field. In the present work, we have used the chiral 
extension |14] (CCDM) of this model to calculate quark matter EOS. The Lagrangian 
density of CCDM is given by 

C(x) = JixW-fd^ - (m + m/ X (x)U 5 ) + (l/2)g^\ a A°(x)}i> 

+f 2 j4Tr(d,Ud»Uj) - l/2mj0 2 (x) 
-(l/4) X 4 (x)(F;,(x)) 2 + (l/2)a 2 v (d, X (x)) 2 - U( X ) (8) 

where U = e tXa<l>a ^ n and U 5 = e lXa<f ' a " ys ^' K , ip(x), A^(x), x(%) and <j)(x) are quark, gluon, 
scalar ( colour dielectric ) and meson fields respectively, m^ and m are the meson and 



quark masses, f n is the pion decay constant, F^(x) is the usual colour electromagnetic 
field tensor, g is the colour coupling constant and A a are the Gell-Mann matrices. The 
flavour symmetry breaking is incorporated in the Lagrangian through the quark mass 
term (mo + m/xUs), with mo = for u and d quarks. So the masses of u, d and s 
quarks are m, m and itiq + m respectively. The self interaction U(x) of the scalar field 
is assumed to be of the form 

aB X 2 (x)[l - 2(1 - 2/a) X (x) + (1 - 3/a) X 2 (x)} (9) 

so that U(x) has an absolute minimum at x = and a secondary minimum at % = 1. 
The interaction of the scalar field with quark and gluon fields is such that quarks and 
gluons can not exist in the region where x — 0- m the limit of vanishing meson mass, 
the Lagrangian of eqn. (8) is invariant under chiral transformations of quark and meson 
fields. 

The calculation of equation of state proceeds as follows. We assume that, in presence 
of nonzero quark/anti-quark densities, the square of meson fields may develop nonzero 
vacuum expectation values (< <fi 2 >) [|TJ], |T8fl . This assumption is analogous to the 



assumption that, in linear a model |15[, the a field acquires nonzero vacuum expectation. 



In ||14|| , it has been shown that this occurs when the quark density exceeds certain critical 
value and one of the effects of nonzero < <p 2 > is that the effective quark masses decrease 
with the increase in < cf 2 >. Thus at large quark densities, one obtains an equation of 
state similar to the equation of state of free quarks and gluons. Further more, we adopt 
mean field approximation to calculate the colour dielectric field (x) in quark matter. 
With these assumptions, the thermodynamic potential for the quark matter is given by, 



-^[k^-mfH^^^-kt]} (10) 



where i — u,d, s. Also, fcj = y//f — m* 2 and this becomes equal to Fermi momentum 
for a s = 0. Here Q is calculated upto second order in quark-gluon interaction. 
In addition, chemical equilibrium under weak decay and charge neutrality imply, 

^s = Hu + He (11) 

and 

2/3n u - l/3n d - l/3n s - n e = (12) 

The baryon density tib = 1/3 J2ii n i) where i — u,d,s and baryon chemical potential is 
defined as hb = Hd + Hu + Hs- 

The mean field values of x an d F$ are calculated by minimising Q. Equivalently, one 
can solve the equations of motion for x, F V1 Fk and F v as obtained from the Lagrangian. 
We find that < n 2 > alone develops non zero value in the quark matter. < K 2 > and 
< rf > remain zero throughout the range of densities considered. This implies that 
strange quark mass remain constant in the medium. On the other hand, u and d quark 
masses change in the medium. The pressure P and energy density E of the quark matter 
are calculated using the relations, P = —Q and E = Q + fiifii. Here, the number 
density of the quark of ith type is given by rii = —(dfl/dfii) 

The parameters of the CCDM are u and d quark mass, strange quark mass, strong 
coupling constant a s , bag pressure B and a. These are obtained by fitting the baryon 
masses. Earlier calculations [[17]] show that these are not determined uniquely by the 



fitting procedure. In particular, it has been found that good fits to baryon masses are 



obtained for 0.6GeV < m GB < 3GeV , m q {u, d) < 125MeV, m q {s) ~ 300MeV and 
£? 1//4 < 150MeV, where glueball mass itigb is defined as vr? GB = 2Ba/a^. The fits are 
better for lower values of mcB, m q (u,d) and B. Since we do not consider scalar field 
excitations, the quark matter EOS does not depend on mas- 

3 Hadron - Quark Phase Transition 

Having determined the equations of state of neutron and quark matter, the phase 
transition point is determined by adopting Gibb's criterion. That is, the point at which 
the free energies ( or pressure ), for given chemical potential, are equal gives the phase 
transition point. In the present context, the baryon and electron chemical potentials ( fis 
and /i e respectively ) are the two independent chemical potentials. However, if we assume 
the charge neutrality, /i e is no more independent. Then the crossing of the pressure curves 



for two phases in P -jj^b plane gives the phase transition point. Glendenning ||11| , on the 
other hand considered a case where, in the mixed phase, neutron and quark matters are 
not charge neutral but the mixture as a whole is. This aspect has been studied further 



by Heiselberg et.al. in ref p0| . We are not considering this situation here. 

For a first order transition, the derivatives of P- [ib curve for the two phases at the 
phase transition point are not equal and the difference in the two derivatives gives the 
discontinuity in the density of the two phases at the transition point. The two phases 
coexist in this range of density. The latent heat of transition is given by the difference 
in the energy densities of the two phases at the critical point. As mentioned earlier, the 
lattice QCD calculations indicate that hadron - quark phase transition may be weakly 
first order or second order. In a calculation, such as ours, we will necessarily have a 
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first order transition since two different models are employed to calculate the properties 
of quark and hadron phases. However, our calculation shows that the phase transition 
can be made weaker by varying the compressibility of the nuclear matter as well as the 
meson- strange baryon coupling within reasonable limits. For example, Table 2 shows 
that the latent heat can be reduced from 154MeV A //m 3 to 96MeV/ fm 3 by changing the 
compressibility from 250MeV to 300MeV. Thus, it is possible to mimick the second 
order phase transition in this manner. 

Our calculation shows that as the nuclear compressibility is increased from 250MeV, 
the chemical potential ( and baryon density ) at which the phase transition occurs 
decreases. ( see Fig. 1) This can be understood as follows. For a given parametrization, 
the nucleon- nucleon repulsive interaction increases with compressibility and hence the 
slope of the P — \xb curve at the transition point i.e. the baryon density decreases. The 
dependence of latent heat on compressibility is given in Table 2. 

As mentioned earlier, the strange quark - meson coupling constant in hadron phase is 
not well known. We have therefore, varied this coupling constant as mentioned in earlier 
section and investigated its effect on the properties of phase transition. We find that as 
the coupling constant is reduced from the universal coupling ( g Ha / 'g Na =l, a = a, oj or 
p and H = A or £~ ) the slope of the P — fiB curve ( Fig. 2 - 5) increase which leads 
to the increase in transition density and as well as reduction in the latent heat (Fig.6 
and Table 2). The EOS becomes softer with decrease in hyperon couplings. This can be 
understood as follows. With decrease in hyperon couplings, it is energetically favourable 
to convert nucleons into hyperons as hyperons do not feel the predominantly repulsive 
force. As a result with decreasing coupling more and more hyperons get populated 
thereby reducing the energy further. 
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As mentioned in the previous section, the quark matter EOS has been obtained 
by using the parameter sets which fit the baryon masses. We find that the transition 
density for hadron quark phase transition decreases rapidly with the decrease in B. 
For I? 1 / 4 < lAOMeV, the transition density is smaller than the nuclear matter density. 
This is unphysical, since it implies that we should have quark matter at nuclear matter 
densities. This happens because the contribution of the dielectric field to the pressure, 
which is negative, is proportional to B. So for smaller B, the pressure of the quark 
matter increase and this makes the quark matter stable at lower chemical potentials 
and densities. This type of behaviour has also been noticed for the bag model EOS 
where one requires B 1 / 4 ~ 150MeV ||. At this point, we would like to note that B can 
not be increased arbitrarily if one insists on a reasonable fit to baryon masses. Thus in 
CCDM, a very restricted parameter sets give reasonable values for transition densities. 
One such set (5 1 / 4 = 152.1MeV, m q{u4) = 91.6MeV, m q{s) = 29A.9MeV a = 36 and 
strong coupling a s = 0.08) has been used in the calculations reported here. 

Finally, let us consider the change in strangeness fraction at the transition point. 
The strangeness fraction in quark and nuclear matter are defined by the ratio strange 
quark density/ baryon density and strange hadron density/baryon density respectively. 
Note that for equal mass quarks and A matter this ratio is unity. We find that generally 
at the transition point the strangeness fraction is larger in quark matter. This is shown 
in Fig. 7, where strange fraction is plotted as a function of baryon density. Note that 
in the coexistence region, the system consist of a mixture of quark and nuclear matter. 
The strangeness fractions in the mixed phase are calculated using the linear relation 
f s (M) = (1 — x)fs{H) + xfs{Q)i where \ is the concentration of the quark matter 
in the mixed phase and f s is the strangeness fraction, M, H, and Q denoting the 
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mixed, hadronic and quark phase respectively. Similar relation is used to calculate the 
corresponding baryon density as well. It is evident from the graph that for a large 
mixed phase region, the jump in the strangeness fraction from hadronic to quark phase 
is larger. 

The presence of considerable amount of strangeness in the hadronic sector in the 
form of hyperons has important consequences on the mechanism of phase transition and 
core of the neutron stars. Some of the possible mechanisms are discussed in ref [ Alcock, 
Farhi and Olinto PT]] ]. The conversion via two flavour has been considered in ref |[22|| , 
where as in ref ||23||, Olinto has started with the assumption that strange matter has 



been seeded into neutron star from outside. On the other hand, presence of considerable 
amount of hyperons near the transition point suggest that inside a cold neutron star 
conversion from neutron to strange quark matter may occur through the clustering of 



A's [^TJ. For a hot neutron star the conversion to strange matter will occur due to 
thermal A's. Also, along with these processes, there will be usual strangeness changing 
weak decay to convert excess d quarks to s quarks. We believe that these points need 
further investigation. 

4 Neutron star models 

Quark cores can exist inside a neutron/hybrid star only for a narrow range of param- 
eter sets in our model. The extent of quark core will be higher for larger compressibility 
and larger hyperon couplings. 

The structure of a neutron star is characterised by its mass and radius. Additional 
parameters of interest are the moment of inertia, the surface red shift z and the rela- 
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tivistic Keplarian rotation period Pk denned as 

P K = 0.026 



(R/km) 3 , . 

r ms] (13) 



\ (M/M Q ) 

as a function of the central density p c of the star. These are important for the dynamics 
and transport properties of pulsars. 

The equations that describe the hydrostatic equilibrium of degenerate stars without 
rotation in general relativity is called Tolman-Oppenheimer-Volkoff (TOV) equations, 
which is given in ref []3], [T^, ^4], 35| . These equations can be numerically integrated, for a 



given central density, to obtain the radius R and the gravitational mass M of the star. 
The moment of inertia I of the rotating neutron star, is also calculated [see ref |3|, [Z5J ] 
for the corresponding central density. To integrate the TOV equations, one needs to 
know the equation of state P(p), for the entire expected density range of neutron star, 
starting from the higher density at the center to the surface density. The composite 
equation of state for the entire neutron star density span, was constructed by joining 
the nonlinear Walecka hadronic equation of state (eqn.6-7) to that given by (a) Negle 
and Vautherin [2B| for the density range 10 14 to 5 x 10 10 gm/cm 3 , (b) Baym, Pethick 



and Sutherland |27) for the range 5 x 10 10 to 10 3 gm/cm 3 and (c) Feynman, Metropolis 
and Teller |2j| for densities less than 10 3 gm/cm 3 . 

The results for star structure parameters are listed in Table 3. Fig.8, and Fig.9 show 
plots of mass vs central density and mass vs moment of inertia respectively. We have 
plotted the curves for hadronic EOS with hyperon couplings (1), (2) and (4) as given in 
section 2.1 and fixed compressibility 300MeV. In case of quark matter, we have used 
the interacting CCDM with the parameter set discussed in section 3. For neutron stars a 
stiffer EOS results in a larger maximum mass and radius 0, 13]. So for heavier stars, one 
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needs larger compression constant K or larger hyperon couplings. Similarly for quark 
stars, one needs larger B or q s . This dependence is reversed for hybrid stars where a 
more repulsive interaction yields lighter stars as can be seen in Fig. 8 and Table 3. This 
is because a stiffer EOS implies lower critical baryon densities in the hadronic sector and 
hence larger quark cores. Since quark EOS yield smaller maximum masses, this reduces 
the maximum mass of the hybrid stars (Fig.8 and Table 3). 

Here we would like to mention that the most accurately determined mass is that 
of PSR1913+ 16 with M/M Q = 1.44 ± 0.003 gj. The observational lower limit on the 
moment of inertia [2~9~| is I = 40M & km 2 ~ 8 x lO M gm/cm 3 . The red shift, which is not 



measured for any neutron star with known mass, seem to lie in the range 0.2 - 0.5. In 
principle, the above observational results on neutron stars should put some constraint 
on the EOS. But in practice, most of the EOS reproduce consistent results as found by 



other authors as well 0, [12|]. This can also be seen from our results for characteristics 
of the maximum mass as given in Table 3. Table 3 contains the results for neutron 
stars models HM J, HM II and HM III corresponding to three hyperon couplings 
((1), (2) and (4), see section 2.1) with compressibility 300MeV and for hybrid stars 
models Hybrid I, Hybrid II and Hybrid III respectively. In our model, the maximum 
gravitational mass for stable non- rotating neutron stars are in the range 1.86M — 1.6M . 
This decrease in maximum mass, is due to softness of EOS. That is because of decrease 
in hyperon couplings as explained earlier. The corresponding radius increases from 
10.69A;m — 11.00km, whereas, red shift and moment of inertia decrease from 0.44 — 
0.32 and 1.76 x 10 45 g cm 2 — 1.52 x 10 45 g cm 2 respectively. Similarly, the hybrid 
star maximum mass corresponding to the three hyperon couplings are in the range 
1.47M — 1.49M . The corresponding variation of radius, red shift and moment of 
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inertia are 12.36km - 12.0km, 0.242 - 0.256 and 1.75 x 10 45 g cm 2 - 1.71 x 10 45 g cm 2 . 
The small variation of hybrid star properties implies that hadronic EOS does not have 
strong influence on the masses and radii of the hybrid star as observed by other authors 
as well 0] . We have also calculated the Keplarian rotation period. The minimal rotation 
periods for neutron stars are in the range 0.66ms — 0.75ms. The corresponding range 
for hybrid stars is 0.93ms — 0.89ms. So the periods for both neutron and hybrid stars, 
in our model, are comparable with the limits obtained by other authors [ffl. Thus, the 
general characteristics of stable star, in our model, are compatible with the observational 
estimates. 

The post glitch data set from the vela pulsar indicates that crust superfluid comprises 
of about 3.4 x 10 -2 of the stars moment of inertia |JIJ. This is lower bound on the 
fractional moment of inertia of the entire crust. We have calculated a = I p /I, where 
I p is the moment of inertia of the pinned superfluid region and I is the total moment 
of inertia of the star. This region is defined as the radial extent corresponding to the 
density 2 x 10 14 — 2 x 10 13 gm/cm 3 . Table 3 shows that the value of a is sensitive 
to the values of hyperon couplings. The purpose of this calculation is to compare an 
observational feature of pulsar glitches with predicted theoretical values of neutron and 
hybrid stars in the present model. 

5 Summary 

We have described the hadron - quark phase transition at zero temperature in the 
frame work of relativistic mean field theory for hadronic sector and chiral Colour Dielec- 
tric model for quark sector. We have studied the hadronic EOS and phase transition by 
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varying the compressibility from 250MeV — 350MeV. Using a simple quark counting 
rule, we get a hyperon coupling which is different compared to those used by other au- 



thors [[TO], FI], . We have given results for two other couplings found in the literature as 
well. We find that critical jjlb and critical ris can be reduced by increasing the hyperon 
coupling or compressibility. Moreover, latent heat can also be reduced by decreasing 
the hyperon coupling or varying the compressibility. This further implies a weaker first 
order transition. In quark sector, we have taken the parameter sets from baryon spec- 
troscopy. We find that only for a very narrow range of parameters in CCDM, one can 
get phase transition. In CCDM, the quark matter EOS depends sensitively on quark 
masses, strong coupling and bag pressure. For _B 1//4 < lAOMeV the transition density is 
less than the nuclear matter density which is unphysical. Here we would like to mention 
that in our calculation, it is not possible to change any of the parameters arbitrarily, as 
it would destroy the fit of baryon masses. 

We find that a considerable amount of strangeness is present in the hadronic phase 
at the transition point. Also there is a large change in the strangeness fraction from 



hadronic to quark phase. Several authors [22] have studied the production rates of 
strangeness during the phase transition from nuclear to quark matter. But as predicted 
by Alcock et.al. [[H]], the rate of transition to quark matter may be much faster, if 



strangeness is already present in the hadronic phase. It will be interesting to study the 
mechanism of transition to quark phase and the rate of strangeness production in such 
circumstances. 

The neutron star characteristics are calculated for different EOS obtained from 
Walecka models and CCDM. It is clear that the mass limits for stars are sensitive to 
the parameter sets. We find that CCDM gives softer EOS compared to bag models and 
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hence lower masses and radii. Hybrid star masses are lower compared to pure neutron 
stars because of the quark cores. We found that the stars properties are not much sen- 
sitive to the composition of stars. Also these characteristics obey the observationally 
inferred limits. So, it is difficult to distinguish between the stars of different composition 
observationally. Hence the most massive pulsars observed so far could possibly either 
a neutron star or a hybrid star. We do not get the layered structure as predicted in 



reference plf . In fact, we find that parameter sets which predict a stable quark phase 
at lower densities, also predict a unstable isospin symmetric nuclear matter at lower 
densities. Here we would like to point out that the nuclear to quark matter phase tran- 
sition can happen through some intermediate phases as well. The nuclear matter may 
go over to a pion or kaon |JI| condensate phase through a second order transition which 
then goes over to quark phase at higher densities. In such case the transition to quark 
phase happens at much higher densities. Some other authors have tried to explore the 



possibility of second order transition |T2j , which we do not consider here. 
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Table 1: Coupling constants for several compresibility K and B/A=-16 MeV, 
p = 0.15 fm~ 3 , a sym = 32 MeV and m* sat /m = 0.8 



K 


(g s /m s ) 2 


(g w /m w ) 2 


(g P /m p ) 2 


b 


c 


(MeV) 


(fm") 


(fm 2 ) 


(fm 2 ) 






250 


9.216 


4.356 


5.025 


0.008209 


0.007385 


300 


8.492 


4.356 


5.025 


0.002084 


0.02780 


350 


7.820 


4.356 


5.025 


-0.004618 


0.05015 



Table 2: Characteristics at the critical point. The columns correspond to compressibility 
(K), hyperon couplings (HC), critical pressure (P c ), critical chemical potential (/i c ), 
energy of the hadronic phase at critical point (E c n), latent heat (L), baryon density 
in the hadronic phase at critical point (nB(cH)) and baryon density width of the mixed 
phase (Aub) 



K 


HC 


Pc 


Mc 


E C H 


L 


nB(cH) 


An B 


(MeV) 




MeV/fm 3 


MeV 


(MeV/fm 3 ) 


MeV/fm 3 


fm~ 3 


fm~ 3 




(1) 


62.2 


1179.5 


490.7 


123.2 


0.47 


0.11 


300. 


(2) 


78.5 


1200.0 


581.2 


96.4 


0.53 


0.09 




(4) 


75.0 


1194.0 


603.7 


58.6 


0.55 


0.06 


250. 


(2) 


269.3 


1428.0 


1236.6 


153.6 


1.0 


0.11 


350. 


(2) 


58.0 


1163.5 


480.8 


128.5 


0.45 


0.11 
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Table 3: Star characteristics for different models. Here IQM for interacting quark 
matter and HM for hadronic matter. 



Pc 


R 


M/M Q 


z 


I 


Pk 


Ip/I 


Model 


(g cm~ 3 ) 


(km) 






(g cm 2 ) 


ms 






3.60xl0 15 


8.14 


1.38 


0.41 


8.22 xlO 44 


0.52 


- 


IQM 


2.50xl0 15 


10.69 


1.86 


0.44 


1.76xl0 45 


0.66 


2.98xl0~ 2 


HM I 


2.50xl0 15 


10.93 


1.72 


0.37 


1.52xl0 45 


0.72 


6.10xl0- 2 


HM II 


2.00xl0 15 


11.00 


1.60 


0.32 


1.52xl0 45 


0.75 


2.68xl0- 2 


HM III 


0.84xl0 15 


12.36 


1.47 


0.24 


1.75xl0 45 


0.93 


2.94xl0- 2 


Hybrid I 
(IQM+HM I) 


0.99xl0 15 


12.94 


1.49 


0.23 


1.72xl0 45 


0.99 


5.4X10" 2 


Hybrid II 
(IQM+HM II) 


1.03xl0 15 


12.00 


1.49 


0.26 


1.71 xlO 45 


0.89 


2.39xl0- 2 


Hybrid III 
(IQM+HM III) 
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FIGURE CAPTIONS 

Figure 1: Pressure vs. chemical potential (a) Interacting quark matter (Q.M.), (b) 
Hadronic matter (H.M.), K= 250MeV , hyperon coupling (HC) (2), (c) H.M., K= 
300MeV, HC (2), (d) H.M., K= 350MeV, HC (2) and (e) Bag model with bag 
constant B 1 / 4 = 152.1MeV 

Figure 2: Pressure vs. chemical potential (a) Q.M. (b) H.M., K= 300MeV, HC (1) 

Figure 3: Pressure vs. chemical potential (a) Q.M. (b) H.M., K= 300MeV, HC (2) 

Figure 4: Pressure vs. chemical potential (a) Q.M. (b) H.M., K= 300MeV, HC (4) 

Figure 5: Pressure vs. chemical potential (a) Q.M. (b) H.M., K= 300MeV, HC= 0.75 

Figure 6: EOS with first order phase transition from H.M. to Q.M. with compressibility 
300MeV (a) HC (1), (b) HC (2) and (c) HC (4) 

Figure 7: Variation of strangeness fraction with baryon density for the EOS in Fig.6, 
(a) HC (1), (b) HC (2) and (c) HC (4) 

Figure 8: Mass vs. central density of neutron (K= 300MeV), quark and hybrid stars 
(a) neutron star (N.S.), HC (1), (b) N.S., HC (2), (c) N.S., HC (4), (d) quark 
star . Corresponding curves for hybrid star are denoted by Hybrid I, Hybrid II and 
Hybrid III respectively. 

Figure 9: Moment of inertia vs. mass of neutron (K= 300MeV), quark and hybrid stars 
(a) N.S., HC (1), (b) N.S., HC (2), (c) N.S., HC (4), (d) quark star . Corresponding 
curves for hybrid star are denoted by Hybrid I, Hybrid II and Hybrid III respectively. 
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